Multilayer network structure enhances the coexistence of competitive species

While the analytical solution for marginal distribution of a stochastic chemical reaction network has been extensively studied, its joint distribution, i.e. high-dimensional master equation, received much less attention. Here we develop novel method computing exact distributions wide class first-order systems in steady-state conditions. The effectiveness our is validated by applying it to four gene expression models biological significance, including with 2A peptides, nascent mRNA, regulation, translational bursting, and alternative splicing.